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EFFECT OF LIMITED AMPLITUDE AND RATE OF FLAP MOTION 
ON VANE -CONTROLLED GUST-ALLEVIATION SYSTEM 

By L. Keith Barker, Daniel J. Crawford, 
and Gene W. Sparrow 
Langley Research Center 

SUMMARY 

An airplane (light transport type) is assumed to be in level flight (no pitching) 
through atmospheric turbulence which has a mean- square vertical gust intensity of 
9.3 (m/sec)^. The power spectral density of the vertical acceleration due to gusts is 
examined with and without a gust-alleviation system in operation. The gust-alleviation 
system consisted of wing flaps that were used in conjunction with a vane mounted ahead 
of the airplane to sense the vertical gust velocity. 

The primary purpose of this study was to examine the change in the effectiveness 
of the gust-alleviation system when the flap motion is limited in amplitude and rate. The 
alleviation system was very effective if no restrictions were placed on flap motion (rate 
and amplitude). Restricting the flap amplitude to 0.5 radian did not appreciably change 
the effectiveness. However, restricting the flap rate did reduce the gust alleviation, and 
restricting the flap rate to 0.25 rad/sec actually caused the alleviation system to increase 
the vertical acceleration above that for the no -alleviation situation. Based upon this anal- 
ysis, rate limiting appears to be rather significant in gust-alleviation systems designed 
for passenger comfort. 


INTRODUCTION 

An airplane flying in turbulent air is sometimes subject to motions which cause 
pilot and passenger discomfort or sickness. For this reason, studies have been made 
of the use of gust- alleviation systems. The basis of one such system involves using a 
vane mounted ahead of the aircraft as a gust sensor, and deflecting various airplane con- 
trol surfaces in response to the vane motion in order to generate forces to counteract 
gust forces. Theoretical studies have indicated that such a system can be very effective 
in reducing the aircraft motion in gusty air. However, such studies generally have placed 
no limit on the amplitude or rate of control-surface motion. (See refs. 1 to 5, for exam- 
ple.) Since practical systems would have some limit on amplitude and rate of motion, it 
is of interest to determine how such alleviation systems are affected by these limitations. 



The present analytical study, therefore, was made to examine the effects of such 
limitations. In order to simplify the analysis, the airplane was restricted to vertical 
translation (no pitching). A signal from the gust- sensing vane is used to move the flaps 
either upward or downward through a linear servomechanism to counteract the effects 
of the gusts. The effectiveness of the gust-alleviation system is evaluated by examining 
the power spectral density of the vertical acceleration of the aircraft. The study was 
restricted to one-dimensional gusts (that is, no spanwise variation of gust intensity) since 
spanwise variations were found to have little effect unless the wing span becomes compa- 
rable with the scale of turbulence (30 meters or more). (See ref. 4.) 

SYMBOLS 

Cg vertical-force coefficient 

3C7 

Cz « = ~da~ ’ per radian 
3C7 

c Z5f ~ “agp P er radian 

E servosystem input signal, rad 

f vertical gust frequency, Hz 

g acceleration due to gravity, m/sec2 

K scalar gain constant in flap servosystem 

K2 gain constant in gust spectral shaping filter 

L scale of turbulence, m 

m mass of airplane, kg 

S wing area, m2 

t time, sec 

V forward-flight velocity, m/sec 

Wg vertical gust velocity (positive upward) , m/ sec 

w 2 mean- square vertical gust velocity, (m/sec) 2 

o 


2 



a 

A a 

°fe 

°V 

6f 

Wl 

Ml 

^NS 

P 

& NS 

r 

'fiwg 
<t >. z 
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Wn 


vertical displacement (positive downward) , m 
angle of attack on wing, rad 

change in angle of attack on wing from trim condition, rad 

angle of attack on wing due to vertical gust velocity, rad 

total angle of attack of vane, rad 

flap deflection, rad 

limit on flap deflection, rad 

limit on flap rate, rad/sec 

damping ratio in flap servo system 

mean estimator of filtered white noise 

air density, kg/m 2 

standard- deviation estimator of filtered white noise 

time required for gust to travel from vane to wing, sec 

one-sided power spectral density of vertical gust velocity, (m/sec2)^/Hz 

one-sided power spectral density of vertical acceleration of aircraft, 
(m/sec 2 ) 2 /Hz 

vertical gust frequency, rad/sec 

natural frequency of flap servosystem, rad/sec 


Primes denote derivatives with respect to distance. Dots over symbols denote 
derivatives with respect to time. Tildes denote functions of distance. 


ANALYSIS 


The airplane motion is one -dimensional (translation in z-direction, no pitching). A 
vane is mounted ahead of the airplane to sense the vertical gust velocity. The vane mea- 
surement is then used to deflect the flaps either upward or downward in order to counter- 
act the gusts and reduce the vertical motions of the airplane. 
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The power spectral density of the vertical acceleration of the airplane is computed 
and is used to examine the effectiveness of the gust-alleviation system. With no restric- 
tions on flap motion, the equations used are linear and a closed-form expression can be 
obtained to relate the power spectral density of the aircraft response to the power spec- 
tral density of the gust. However, when velocity and/or displacement limits are speci- 
fied for flap motion, the equations become nonlinear and the power spectral density of 
the aircraft response must be estimated from the computed time history of the aircraft 
response to the random gust. Such time histories were generated by an analog computer 
and, subsequently, analyzed by a digital spectral analysis procedure. 

Equation of Motion 

The equation of motion used in this study for an aircraft free to move only in the 
vertical direction is given by 

mz(t) = £^[Cz a AQ'(t) + 5f(t)j (1) 

where Aff(t) is defined by 

A«(t) = ^ + «g(t - r) = ^ + Wg y~ T) (2) 

Here Wg(t - r) is the vertical gust velocity which hit the vane r seconds ago 
and is presently acting on the wing. Hence, equation (1) can be written as 

mz(t) =^y<|cZaj^ + Wg( y + C Z6f 6f(tj> (3) 

This equation can also be written as 

z(t) + ajz(t) = -ajWgft - t) + a2 6 f( t ) ( 4 ) 

where 

a l = ~P §5T C Za (5) 

and 

a 2 ’ >> 1st Cz «t (6> 

It is assumed that the flap servomechanism can be described by a second-order 
differential equation. (See ref. 2.) 
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6f(t) + 2?« n 6f(t) + w n 2 6f(t) = -Ka) n 2 E(t) 


(7) 


Two reasonable servosystem inputs are 


E(t) = a v (t) = 


wg(t) 

V 



( 8 ) 


and 

E(t) = Og(t) = = o^t) - ^ (9) 

Use of equation (8) means that the flap is moved in response to the total deflection of the 
vane from its trim position. Use of equation (9) means that the flap moves in response 
to the true gust angle of attack (that is, the vane angle corrected for airplane motion). 

The latter case would be more difficult to mechanize than the former case. A prelimi- 
nary analysis was made to examine each of the two types of control inputs. (See appen- 
dix A.) The results showed that rate limiting of the flap motion had similar effects on 
the power spectral density of the response of the aircraft to random gusts. Somewhat 
better alleviation was obtained by using just the gust angle of attack (eq. (8)) for the servo 
input. This result appears to be caused primarily because of a limit-cycle oscillation of 
the flap which is induced by the rate limits with the z/V feedback of equation (8). 

Since the effects of rate limiting appeared to be about the same for both types of 
servo inputs, and to avoid any possible problems with limit cycling, the remainder of the 
study was restricted to use of equation (9) for the servo input. 

When equation (9) is substituted into equation (7), the result is 

6f(t) + 2£o> n 5 f (t) + co n 2 6f(t) = -Ka) n 2 (10) 


The gain constant K used in equation (10) was selected to be 


C Z6f 


(ID 


This value corresponds to the optimum value of K when the flaps move in phase with a 
vane that is mounted on the wing aerodynamic center (r = 0), and there are no constraints 
on the flap motion. 


If the time variation of the vertical gust velocity Wg(t) is specified, equations (3) 
or (4) and (10) can be used to compute a time history of the vertical acceleration of the 
aircraft. From this time history, the power spectral density of the vertical acceleration 
can be computed. 
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Generation of Vertical Gust Velocity 

As stated previously, placing limits on the flap-motion velocity and amplitude 
resulted in a set of nonlinear equations of motion for the flap and aircraft response and 
in turn, meant that the power spectral density of the aircraft motion had to be determined 
from the computed time histories. 

To be representative of atmospheric turbulence, it was required that the time his- 
tory of the vertical gust velocity be normal with a zero mean. The following power spec- 
tral density was assumed: 


where 



4a 3 w g 2 
1 + (a 3 o)) 2 


a 3 


L 

2V 


( 12 ) 


(13) 


This power spectral density is based on an approximation of the vertical gust correlation 
function and closely represents turbulence measurements made in wind tunnels. (See 
ref. 3.) 

The techniques for obtaining turbulence having this power spectral density, for 
including the time lag r of equation (4), and for obtaining the power spectral density 
of the aircraft response are given in appendix B. 


Aircraft Characteristics 

The aircraft mass, aerodynamics, and flight conditions used in this study corre- 
spond to those of the light transport airplane of reference 2 and are as follows: 


m, kg 3629 

S, m 2 32.42 

CZq;) Per radian -5.30 

p (sea level), kg/m 3 1.223 

Cz 5 j, Per radian -0.80 

V, m/sec 67 

The two turbulence parameters used in this study are: 

L, m 610 

w g 2 , (m/sec) 2 9.3 
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The present study covered the range of parameters shown in the following table: 


Flap -motion limits 

w n , rad/ sec 



|(5f) L |, rad 

|(fif) L |, rad/sec 


oo 

OO 

10 

0.7 

0.50 

oo 

10 



.50 

0.25 

4, 10 



1 

.50 

1 



l 

1.00 

i 




RESULTS AND DISCUSSION 


Most of the results of this study will be in the form of the power spectral density 
of z, that is, (p£, as a function of frequency f for different flap-amplitude limits 

|(6f) L | and different flap-rate limits |(6f) L |. The power spectral density curves pre- 
sented in this study can be scaled for a wider range of flight conditions as shown in 
appendix C. 


Basic Configuration With No Alleviation 

The simulation power spectral density of the aircraft vertical acceleration is shown 
as a function of gust frequency in figure 1. This curve differs slightly from the corre- 
sponding theoretical curve. Since the power spectral density curves for the nonlinear 
cases must be obtained by the simulation, the simulation curve of figure 1 is used as the 
basis against which gust alleviation will be measured. 

Effectiveness of Flap With No Limits on Amplitude and Rate of Deflection 

The effectiveness of the flap with no limit placed on the amplitude or rate of motion 
is shown in figure 2. With the gust sensor located at the wing aerodynamic center (Vr = 0) 
or at VY = 5.22 meters ahead of the wing aerodynamic center, the flap can eliminate 
most of the acceleration due to gusts. The alleviation is slightly better with the sensor 
ahead of the aerodynamic center since it provides a time difference that allows the flap 
to move toward its commanded position as the gust moves from the sensor to the flap. 

Effectiveness of Flap With Limits on Motion Amplitude 

The change in flap effectiveness caused by limiting the amplitude of flap travel is 
shown in figure 3. It was believed that an amplitude limit of | (®f ) l | = 0. 5 rad would be 
practical for most transport type airplanes. Placing this limit on flap amplitude did not 
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reduce the flap effectiveness in alleviating the gust appreciably for either of the two sen- 
sor positions. (Compare figs. 2 and 3.) 


Effectiveness of Flap With Limits on Both Amplitude and Rate of Deflection 


The change in flap effectiveness caused by limiting the rate of deflection of 
the flap is shown in figures 4 and 5. The results show that a limiting-rate value of 
(6 f ) L | = 1 rad/sec does not affect the ability of the flap to provide good gust allevia- 
tion. Reducing the rate limit to lower values diminishes the alleviation. With a rate 
limit of | (^)l | = 0-25 rad/sec, the flap actually reinforces the gust (improper phasing), 
and the power spectral density of the aircraft acceleration is worse than that for the 
no-alleviation situation. Reducing the natural frequency of the flap system from 
w n = 10 rad/sec to ca n = 4 rad/sec (compare figs. 4 and 5) reduces the maximum 
gust alleviation which can be obtained with a rate limit of |(6f)jJ = 1 rad/sec. 


The effects of rate limiting are summarized in figure 6. This figure shows the 
percent reduction in the root mean square of the vertical acceleration z as a function 
of rate limiting 




, natural frequency co n , and sensor position Vr. The effects of 
the various parameters on the root mean square of z have the same trends as they did 
on the power spectral density of z. 


Based upon this analysis, rate limiting appears to be rather significant in gust- 
alleviation systems designed for passenger comfort (lower frequency range); however, 
because of the limited nature of the study, further analysis is warranted. 


CONCLUDING REMARKS 


An airplane (light transport type) is assumed to be in level flight (no pitching) 
through atmospheric turbulence which has a mean square vertical gust intensity of 
9.3 (m/sec)2. The power spectral density of the vertical acceleration due to gusts is 
examined with and without a gust-alleviation system in operation. The gust-alleviation 
system consisted of using wing flaps in conjunction with a vane mounted ahead of the air- 
plane to sense the vertical gust velocity. 

The primary purpose of this study was to examine the change in the effectiveness 
of the gust-alleviation system when the flap motion is limited in amplitude and rate. The 
alleviation system was very effective if no restrictions were placed on flap motion (rate 
and amplitude). Restricting the flap amplitude to 0.5 radian did not appreciably change 
the effectiveness. However, restricting the flap rate did reduce the gust alleviation, and 
restricting the flap rate to 0.25 rad/sec actually caused the alleviation system to increase 
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the vertical acceleration above that for the no-alleviation situation. Based upon this 
analysis, rate limiting appears to be rather significant in gust-alleviation systems 
designed for passenger comfort. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., February 18, 1972. 
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APPENDIX A 


FLAP SERVOSYSTEM INPUT E(t) 

Two expressions were examined for the flap servosystem input; namely, 

E(t) = ag( t) = ^ = a v (t) - M (Al) 

and 

E(t) = «v(t) = + (A2) 

Equation (Al) represents the angle of attack on the vane due to the gust; equation (A2) 
represents the total angle of attack measured by the vane. 

For a random time variation of the vertical gust velocity Wg(t), equations (3) or (4) 
and (7) can be used in conjunction with either equation (Al) or (A2) to generate a time his- 
tory of the vertical acceleration of the airplane. From this time history, the power spec- 
tral density of the vertical acceleration can be computed. The reason for computing the 
power spectral density directly from the time history is that limits placed on the flap 
motion cause the equations to become nonlinear so that the usual closed-form expression 
cannot be used. 

The power spectral density of the vertical acceleration is shown in figure 7 for 
the two types of input given by equations (Al) and (A2). The natural frequency of the flap 
servo system was w n = 4 rad/sec, and the vane location was Vr = 5.22 meters ahead of 
the wing aerodynamic center. The flap amplitude limit was 0.5 radian. The flap rate 
limits used were 1 rad/sec, 0.5 rad/sec, or 0.25 rad/sec. 

From figure 7 it appears that better gust alleviation is attained by using equation (Al) 
rather than equation (A2), since the power spectral density curve associated with equa- 
tion (Al) is generally below the curve associated with equation (A2). As the flap rate limit 
is made more stringent (lowered), this difference becomes more pronounced. The reason 
the power spectral density associated with equation (A2) is amplified appears to result 
from a limit-cycle oscillation of the flaps which is induced by the rate limits with the 
z/V feedback required in using equation (A2). (See fig. 8.) 
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APPENDIX A - Concluded 


The following table shows the root -mean -square values of z corresponding to the 
various power spectral density curves: 


Flap rate limit, 
|(Sf) L |> rad/sec 

Flap amplitude limit, 

Wl|> rad 

Root mean square of z 
based on - 

Equation (A2) 

Equation (Al) 

0.25 

0.5 

2.1 

2.0 

.50 

1 

1.7 

1.4 

1.00 

l 

1.4 

1.3 


As expected, the values associated with equation (A2) are larger than those associated 
with equation (Al). 

Although the results appear to be better when the input given by equation (Al) is 
used, it would appear that equation (A2) would be easier to mechanize. The limit-cycle 
oscillation might possibly be avoided by using some type of nonlinear flap control law. 
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APPENDIX B 


COMPUTER SIMULATION AND SPECTRAL ANALYSIS PROGRAM 

A block diagram of the gust-alleviation system is shown in figure 9. The vertical 
gust velocity Wg(t) is sensed by the vane. Later (after transport delay), this same gust 
velocity Wg(t - r) hits the wing. The flap, meanwhile, moves in response to the gust 
presently sensed by the vane Wg(t). 

The mathematical models of the flap servo (eq. (10)) and the aircraft (eq. (4)) were 
implemented on a general-purpose analog computer which was operated in conjunction 
with a digital computer operating in the real-time mode. The functions of the analog 
computer were to generate the gust disturbance velocities and to solve the equations of 
motion of the flap and aircraft. The functions of the digital computer were to provide the 
transport delay which appears in equation (4) and to calculate the power spectral density 
of the gust input and aircraft response. 

Atmospheric Turbulence Generator 

Figure 10 is a block diagram of the vertical gust generator. It consists of a noise 
source, a spectral shaping filter, and a circuit to insure that Wg has a mean of zero and 
the desired variance Wg2. The noise source was a variable -frequency, constant- variance 
device which gave a continuous Gaussian signal with a mean-square voltage of 10, with a 
selected half -power point at 15 Hz. The spectral-shaping filter was a first-order filter 
with half-power point at about 0.03 Hz (that is, V/ wL). Thus, the spectrum of the noise 
source, with the half-power point at 15 Hz, was much flatter than the spectrum after the 
filter and, consequently, had an almost white spectrum with respect to the filter. The 
static gain of the filter K2 was somewhat arbitrary because an automatic gain circuit 
controlled the mean and variance from the desired simulated gust. 

Figure 11 shows the method of estimating the mean and standard deviation of the 
filtered gust signal. The integrators were allowed to operate at least 1000 seconds 
before being used in the simulation. The purpose was to achieve good statistical esti- 
mates of the mean and standard deviation prior to any analysis of the aircraft response. 

The integrators were not reset for several simulation runs. This estimation technique 
required the generation of the function t - * for t ^ 1000 seconds. This function was 
implemented with the circuit shown in figure 12. 

Aircraft Dynamics 

The diagram for simulation of the flap motion is shown in figure 13. Either of the 
two limits (rate or amplitude) could be selected. The limiter on flap position drives the 
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APPENDIX B — Continued 


flap rate to zero when the flap is at its limits of travel; a conventional velocity limit was 
used for flap rate. The two limiters acted independently. 

The variable time delay was implemented by use of a potentiometer, two analog-to- 
digital converters, a digital-to-analog converter, and a data buffer in the digital computer. 
Although the potentiometer could be varied continuously, the corresponding time delay 
varied in discrete steps of 1.0 sample period (0.03125 second) from a value of 1.5 sample 
periods (0.046875 second) to a value of 32.5 sample periods (1.015625 seconds). The digi- 
tal computer was bypassed when simulating the zero-transport-delay condition. The 
dynamics of the aircraft were simulated as shown in figure 14. 


Spectral Density Calculations 

The power spectral densities of two quantities were desired: that of the aircraft 
acceleration (for evaluation of the gust-alleviation system) and that of the gust input (to 
ascertain that the gust characteristics were shaped as desired). Both quantities were 
computed in the same manner. The following description of the mechanization of the 
computer for power spectral density computations is written in general terms. 

The one-sided spectrum of a stationary random process x(t) is given by 

0 X ( w) = 2 f 0 X ( T ) e~ i0JT dr (Bl) 


where 0 x (t) is the autocorrelation function of x(t). IE x(t) is also ergodic, the auto- 
correlation function of x is defined by 


9 x (t) = 


T 

lim i f x(t) x(t - t) dt 
T—°o 1 J 0 


(B2) 


To insure that the integral in equation (Bl) is bounded, the auto covariance function 
i// x (t) is used instead of the autocorrelation function 0 x (t). If x(t) has a zero mean, 
these two functions are equivalent. The autocovariance function is defined as 


X (t) = 0 x (t) - x2 (B3) 

where x^ is the mean- square value of x(t). In the present study the autocovariance 
function was approximated by the summation 

N 

^ x (t) = * , Y x ( k At) x (k At - r) - i 2 (B4) 

JN - l l—/ 

k =1 

where x is the estimate of the mean value of x, and t - l At. In the summation, 

N + 1 represents the number of sample points on the time history. 


I 
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APPENDIX B - Continued 


Since the autocorrelation 0 x (o>) is an even function, 

r°° 

<p x (oo) = 4 \ I^ x (t) cos (jot dt 
J 0 

In the present study 0„(w) was approximated by the summation 


(B5) 


where 


$ x ( co) = 2 At 


Atp 


P-1 


i// x (0) + 2 ^ x (j At) cos(j Aten) + i// x (p At) cos(p Aten) 

j=i 


(B6) 


and p + 1 is the number of points on the power spectral density curve, and the index k 
is called the harmonic number. This curve was smoothed by using a Hanning filter. (See 
ref. 6.) 

The sample period At used in this study was 0.125 second. Approximately 21 min- 
utes of running time T r was required for each simulated run. This time was determined 
by the desired variance on the statistical estimates and also by the resolution of the spec- 
tral density curves. The criterion used to select the running time is 

T r = -i- (B7) 

/V 2 

where 


/3 = _ 1 - 

p e p At 

and e 2 is a measure of the variation of the estimate and is defined by 
_ var[0 x (co)] 


<p x 2 (co) 


(B8) 


(B9) 


In the present study, the values used for /3 and e were 0.08 and 0.1, respectively. 

When no limits are placed on 5f and 6f, a closed-form expression can be obtained 
for the power spectral density of the aircraft acceleration. The expression is 


= 


(aiw) 2 


cn 2 + a^ 2 

jo> 4 + 2w n 2 (2£ 2 - l)cn 2 


- 2(cn n 2 - cn 2 )cos rcn 


- 4£cn n cn sin TO) +l)0 W(r (o>) 
where aj is defined by equation (5). 


(BIO) 
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APPENDIX B - Concluded 


In order to check the accuracy of the computer simulation and the power spectral 
density computations of the computer, results for <^(a>) of equation (BIO) were com- 
pared with those of the computer. The results are shown in figure 15 and show good 
agreement. 
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APPENDIX C 


APPLICATION OF RESULTS TO WIDER RANGE OF FLIGHT CONDITIONS 

The purpose of the following discussion is to present a procedure for interpreting 
the power spectral density curves of this study for different forward -flight speeds and 
gust intensities. 


Scaling Results for Changes in Forward- Flight Speed 

Taylor's hypothesis which assumes the turbulence to be "frozen" in space is often 
employed in studying the effects of rough air on airplanes. In this manner, the vertical 
gust velocity is essentially a function of position in space. Let A be a space parameter 
so that 


A = Vt 


(Cl) 


and 


Equations (3) and (10), written in terms of the distance A, become 

+ c Z 5f 5f(^ 

Wp-(A) 


mz"(A) = ^<Jcz a fz'(A) + Wg(A y V — 


(C2) 


6f”(A) + 2?S2 m 5 f '(A) + fi m 2 6 f (A) = -K £l m 2 


V 


(C3) 


The primes denote differentiation with respect to A, the tilde signs denote functions of A, 
and 

(C4) 

The limits on the flap amplitude motion in space and time are the same. However, 
the limit on flap rate must be scaled by the velocity V as 

(V)l-^Wl < C5 > 

The power spectral density of the space function z"(A) is related to the power 
spectral density of the time function z(t) by 


<%"( A) (fi) = ^3 ^z(t) (a)) 


(C6) 
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APPENDIX C - Continued 


where 

n = | (C7) 

Equations (C6) and (C7) are derived in appendix D. 

Equations (C2) and (C3) show that if the gust intensity is scaled in space in direct 
proportion to velocity, then the vertical motion will remain the same when plotted against 
X, provided Vt, and co n /V are held constant. This statement also means that the 
power spectral density remains the same when plotted as a function of S2. 

Consider two forward-flight speeds V and Vj = CV (V is the flight speed used 
in this study). Figure 16 shows the various relationships which are involved in going 
from to 0~„ ( X )M and then going from to i( t )( w l)- The sub- 

script 1 denotes motion corresponding to Vj. 

The relations shown in figure 16 can be used to obtain the desired transformation 
equations between and l(t)( w l)’ These equations are 


/ \3 

0z,l(t)H) =(^r) 0z< w ) = c3< M w ) 

(C8) 

(&f,l) L = vH^)l = C (®f)L 

(C9) 

Vi _ 

COj = CO = Cco 

(CIO) 


In brief, assume the gust angle of attack to be a constant space function. Let V change 
by a constant factor C. If u> n/ /V, and Vt are held constant and changes by 

the factor C, then the power spectral density curves of this study apply by simply chang- 
ing the ( GtJ ) an< * 00 ( or ^ scales by C^ and C, respectively. 


Scaling Results for Changes in Gust Intensity 

The power spectral density curves of this study can also be scaled for different gust 
intensities while the same forward speed V is kept. Suppose the vertical gust intensity 
is scaled by a constant factor B. Let the resulting motion be denoted by the subscript 1. 
Equations (3) and (10) then can be written as 



(Cll) 


(C12) 
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Assume that the flap amplitude varies directly by the same factor B as the gust 
intensity. Thus, 

6f }1 (t) = B6f(t) (C13) 


and 


6f,l(t) = B6f(t) (Cl 4) 

By comparing the forms of equations (Cll) and (C12) with the forms of equations (3) and 
(10) under this assumption, it can be seen that 

Bz(t) = zj(t) (Cl 5) 

From equation (Cl 5), it then follows that 

= B 2 (^(w) (C16) 

The limits on the flap motion follow from equations (C13) and (C14) as 

Ml = B («f) L < C17 > 

and 

(5f,l) L = B (^)l (C18) 


The desired transformation equations are equations (C16), (C17), and (C18). 
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APPENDIX D 


DEVELOPMENT OF THEOREM AND APPLICATIONS 
Theorem Development 

The following theorem, developed in the course of this study, is very useful in 
relating power spectral density functions when certain conditions exist. 

Theorem: Consider two stationary random processes f(x) and g(y), where 


y = C 1X 


(Dl) 


and Cj is a constant. If the processes are ergodic and 


f(x) = C 2 g(y) 


where C 2 is a constant, then it can be shown that 


(D2) 


^f(x) (Ax ) = C 2 2 ^ g (y)( A y) 


(D3) 


where i// f ^(Ax) and i// g ^(Ay) are the autocorrelation functions of f(x) and g(y), 
respectively. Also, 


, C ? 2 

%(x) K) = 0g(y) ("y) 


provided the frequencies u> x and Wy are related by 

.. _f^x 
y " C 1 


(D4) 


(D5) 


where ( w x) ancl ^g(y)( w y) are the power spectral density functions of f(x) and 

g(y). 

Proof: The autocorrelation function for f(x) is defined by 

^f(x)( Ax ) = lim ^ J f ( x ) f(x + Ax) dx (D6) 

Since f(x) = C 2 g(y) and Ay = Cj Ax, 


^ f( x)( Ax ) = s(y) g(y + ^y) d y = c 2 2 ^ g ( y )( A y) 


(D7) 
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APPENDIX D — Continued 


IIIlHIIIIlllllllllil 


The power spectral density of f(x) is defined by 

= X" ^£(x)( ix > e(-‘“ xAx) d(Ax) 
Substituting equation (D7) into equation (D8) gives 

K) - c 2V g(y) (Ay) d(ix) 

Since Ay = Cj Ax, equation (D9) can be written as 

W Wx > = ^ £„ W Ay) e ( " IWy4y) d < A y> 

where 


By definition, 

■>g(y)(“y) = L ^gCy) (Ay) d( " 1 “ y Ay) “(Ay) 

From equation (D12), it can be seen that equation (DIO) can be written as 

q 2 

^f(x)K) = “cT 0 g(y)^y) 


(D8) 


(D9) 


(DIO) 


(Dll) 


(D12) 


(D13) 


Applications 


As an example in the use of the theorem, let z(t) = V^z"(x), where X = Vt. In this 
case, f(x)=z(t), g(y) = z"(X), C 2 = V 2 , and Ci = V. Let u> x = w and w y = S2. By 
direct substitutions into equations (D4) and (D5), it follows that 


^z(t) (w) =v3 ^z”(x) (J2) 


where 



Another interesting example is the relationship between the power spectral densities 
of the vertical gust velocity written as a function of space X and time t. Here 

w g (t) = w g (X) 
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where 


X = Vt 

In this case, f(x) = w g (t), g(y) = w g (X), Cj = V, and C 2 = 1. Let u> x = a) and 
u)y = S2. Then, application of the theorem shows that 

*w g (t)("> - 7 *s g w< n > 

where 


It should be noted that some authors define the power spectral density function in a 
slightly different manner than that indicated by equation (D8) or (D12). Sometimes a fac- 
tor of 1/tt or 1 /2 7T will precede the integral. The results of the theorem, however, 
still apply in each case. 
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Power spectral density, </>.., (m/sec 2 ) /Hz 



Gust frequency, f, Hz 

Figure 1.- Power spectral density of the vertical acceleration 
response to v ertic al gusts for the basic airplane (without 
alleviation). w„^ = 9.3 (m/sec)^; l = 610 m. 




Power spectral density, 



(a) Gust sensor located at wing aerodynamic center. Vr = 0. 

Figure 2.- Effectiveness of alleviation system with no limits on flap motion 

o) n = 10 rad/sec. 



Without alleviation 
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Gust frequency, f, Hz 

(b) Gust sensor located at Vt = 5.22 meters. 
Figure 2.- Concluded. 
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(b) Gust sensor located at Vr = 5.22 meters. 
Figure 3.- Concluded. 







(a) Gust sensor located at wing aerodynamic center. VY = 0. 


Figure 5.- Effectiveness of alleviation system 
flap motion. w n = 4 rad/sec; | 


with limits on rate 
| = 0. 50 rad. 


of 
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Flap rate limit, l(Sf) L l, rad/sec 



Figure 6.- Effect of flap rate limit on root mean square of vertical acceleration. 

|(«)l| = °' 50 - 
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Power spectr 



Figure 7.- Power spectral density of vertical acceleration resulting from the use of 
equations (Al) and (A2). u> n = 4 rad/sec; Vr = 5.22 meters; = 0.5 rad. 
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Figure 10.- Turbulence generator. 
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Figure 11.- Estimation of mean and standard deviation. 
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Figure 12.- Generation of (1/t) for t g 1000 seconds. 
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Figure 15.- Comparison of analytical and simulation solutions (no limits). w n = 4 rad/sec; £ = 0.7; r = 0.078 sec. 




Figure 16.- Relationships between power spectral densities when gust intensity varies 
with velocity and Vr, o> n /V, an< * (^L are ^eld constant. 
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